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ABSTRACT 

We have mapped the northern area (30' x 20') of a local group spiral galaxy M33 in ^^CO(J= 1-0) 
line with the 45-m telescope at the Nobeyama Radio Observatory. Along with Ha and Spitzer 24- 
/zm data, we have investigated the relationship between the surface density of molecular gas mass and 
that of star formation rate (SFR) in an external galaxy (Kennicutt-Schmidt law) with the highest 
spatial resolution ('^80pc) to date, which is comparable to scales of giant molecular clouds (GMCs). 
At positions where CO is significantly detected, the SFR surface density exhibits a wide range of over 
four orders of magnitude, from Ssfr ^ 10^^" to ^ 10"^ Mq yr~^ pc~^, whereas the values are 
mostly within 10-40 M© pc~^. The surface density of gas and that of SFR correlate well at a ~l-kpc 
resolution, but the correlation becomes looser with higher resolution and breaks down at CMC scales. 
The scatter of the Esfr~Sh2 relationship in the ^ 80-pc resolution results from the variety of star 
forming activity among GMCs, which is attributed to the various evolutionary stages of GMCs and 
to the drift of young clusters from their parent GMCs. This result shows that the Kennicutt-Schmidt 
law is valid only in scales larger than that of GMCs, when we average the spatial offset between GMCs 
and star forming regions, and their various evolutionary stages. 



1. INTRODUCTION 

Understanding the formation process of massive stars 
in galaxies is one of the key issues in modern astronomy 
because massive stars play essential roles in the evolution 
of galaxies. Because stars are formed from molecular gas, 
the study of molecular gas in galaxies provides us with 
important clues to the fundamental physical processes of 
massive star formation in galaxies. Observational studies 
of galaxies on global scales have shown that the surface 
density of star formation rate (SFR) and that of cold gas 
obey a simple relation, Ssfr k S^^^g, where Ssfr is the 
SFR per unit area and Egas is the surface density of gas. 
This relationship is known as the Kennicutt-Schmidt law 
(hereafter, K-S law: ISchmidtl 119591 : lKennicutt][l998[ ). It 
has been shown that the disk-averaged surface density of 
SFR is much better correlated with that of molecular gas 
Sh-^ than with that of total (H1-I-H2) gas Shi+Hs (e-g-, 
I Wong fc Bliti[20Ql . 
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Because almost all of the existing studies on the K- 
S law have been carried out based on CO observa- 
tions of kpc-scale r esolution or disk-averaged data (e.g., 
iKomugi et al.l[2005l ). not much is known about the valid- 
ity of the K-S law for smaller molecular structures such 
as giant molecular associations ( GMAs) and giant molec- 
ular clouds (GMCs). Recently, iKennicutt et al.l ()2007f ) 
investigated K-S law in M51 down to a linear scale of 
500 pc and cloud mass scales of IO^-IO'^Mq. They found 
that the n onlinear K-S law can be extended down to 
that scale. iBigiel et"al1 ()2008[ ) revealed that the K-S law 
is applicable at the 750-pc scale for seven spiral galax- 
ies, and that it hol ds down to the 250-pc scale in M51. 
iVerlev et al.l ()2010l ) have shown that a loose correlation 
exists even in the 180-pc scale in M33, wh ere a strong 
corre lation was found on the global scale (jHever et al.l 
|2004| ). Because most of the molecular gas is confined 
within molecular clouds and virtually all of the GMCs 
are sites of star formation in the Milky Way, GMCs play 
an important role in star formation. It is important to 
understand how star formation at GMC scales is linked 
to the K-S law, which is valid for kpc-scales. 

One way to address this issue is to conduct a high- 
spatial-resolution mapping of entire molecular gas disks 
in the nearest galaxies. The recent improvements in the 
resolution, sensitivity, and observation efficiency of radio 
telescopes have made such observational studies feasible. 

M33 is the best target for this purpose . It is one of the 
neare st spiral galaxies (D = 840 kpc; iFreedman et al.l 
[T99I in which individual GMCs can be resolved using 
present-day instruments. Furthermore, because its disk 
is moderately face-on (i = 51°), it is easy to study the 
molecular gas in relation to star formation. 

Because of these desirable characteristics, M33 has 
been the subject of many large-scale studies of GMCs. 
Pioneering surveys of the central 1.7-kpc region of M33 
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were conducted with the 12-m t elescope at the Nationa l 
Radio Astronomy Observatory (jWilson fc Scovillill989D 
and the Owens Valley Radio Observatory interferometer 
(jWilson fc Scov illc 1990). They revealed that the prop- 
erties of GMCs (velocity width, diameter, peak bright- 
ness temperature, and mass) were very similar to those 
in the Milky Way Galaxy. More recently, a CO ( J= 1-0) 
all-disk survey was performed w ith the Berkeley-Ihinois - 
Maryland Association (BIMA) (jEngargiola et al.l [20031 ) . 
which identified 148 GMCs with masses over 1.5 x 10^ Mq 
across the star-forming disk. The BIMA survey recov- 
ered only 20% of the CO flux in M33 be cause of the 
missi ng short spacings of the int e rferom eter (|Hever et al.l 
120041 ): thus, iRosolowskv et al.l (|2007| ) performed addi- 
tional CO ( J= 1-0) observations with the 45-m tele- 
scope at the Nobeyama Radio Observatory (NRO), and 
combined them with the existing data obtained from 
BIMA and the 14-m telescope at the Five College Ra- 
dio Astronomy Observatory (FCRAO). They obtained 
the CO (J= 1-0) map in the inner ^15' x 15' (-3.7 kpc) 
region at a resolution of 20" (^^80 pc). 

However, these observations are still insufficient to un- 
derstand the overall nature of the GMCs in M33 because 
they do not cover the entire disk, and they possess in- 
homogeneities due to multi-field interferometric observa- 
tions. Accordingly, a homogeneous wide-area single-dish 
mapping with a high spatial resolution and high sensi- 
tivity is required in order to determine the molecular gas 
properties of the entire disk. In order to achieve this 
purpose, it is effective to adopt the on-the-fly (OTP) 
mapping technique using a single-dish telescope hav- 
ing a large aperture. We performed an M33 all-disk 
^^CO(J=1-0) survey with the 45-m telescope, as an 
NRO legacy project to study the properties and evolution 
of GMCs and star formation within the entire galaxy. In 
this paper, we investigate the K-S law in the GMC scale, 
which is the highest resolution (~ 80 pc) to date, as an 
initial result of this project. 

2. OBSERVATIONS AND DATA 

We observed the ^^CO {J— 1-0) lines toward the north 
in a 30' x 20' (~ 7.3 kpc x 4.9 kpc) region of M33; the 
observation region covered the galactic center and the 
three major Hii regions: NGC 604, NGC 595, and IC 133. 
The mapping area is shown in Fig. 1. 

The observations were carried out from January to 
April, 2008, with the 45-m telescope at NRO. The beam 
size of the telescope (HPBW) was 15" for the rest fre- 
quency of the 12CO(J=1-0) line (115.271204 GHz). 
We used the 25-BE am Array Receiver System (BEARS: 
iSunada et aL|[2000l ). Because BEARS is a double side- 
band receiver, we measured the sideband ratios of each 
beam of BEARS by observing a bright standard source 
NGC 7538 using both BEARS and a single-beam re- 
ceiver SI 00 equipped with a single sideband filter. The 
errors in the scaling factors of BEARS were smaller than 
~20%. The main beam efficiency measured with SlOO 
was 77MB = 0.32 ± 0.02 at 115 GHz. The typical system 
temperature was ^ 700K (SSB). We used digital autocor- 
relators as backends. The spectrometers have 1024 chan- 
nels covering a bandwidth of 512 MHz with a frequency 
resolution of 1.0 MHz, which corresponds to 1332 km s~^ 
and 2.6kms~^. The chopper- wheel technique was em- 
ployed to calibrate the antenna temperature T^. Here- 



after, all the CO intensity measurements are mentioned 
in the Tmb (=?a/??mb) scale. 

The observatio ns were performed using the OTF ob- 
servation mode (jSawada et al.l 120081 ). Throughout the 
observations, the data were sampled every 0.1 s. The 
emission-free "OFF" points were taken at an offset of 
30' from the center. We checked the pointing accuracy 
every hour with a 5-point observation of a SiO maser 
source IRC-F30021 with a 43 GHz SIS receiver (S40). 
We excluded all the data in the observation sequence 
that had pointing errors larger than 7". 5 or wind speed 
V > \Q ms~^ to avoid systematic intensity loss due to 
pointing errors. After these data screenings, the net in- 
tegration time was determined to be approximately 130 
h. 

The data reductions were made with the OTF reduc- 
tion software package NOSTAR, which was implemented 
by the NRO. The data were convolved with a Gaussian- 
tapered Jinc function to create the spectral data cube. 
The grid spacing was taken as 7". 5, and the spatial res- 
olution of the final map was 19". 3. Finally, the baselines 
were fitted and subtracted from the data cube. The ve- 
locity resolution of the final cube is 2.5 kms~^. The rms 
noise in a velocity channel was 130 mK. The rms noise 
in the integrated intensity map was 1.6Kkms"^, which 
corresponds to a mass sensitivity of 4.8 Mq pc~^. 

3. RESULTS 
3.1. CO (J =1-0) Map 

Fig. 1 shows the CO(J=1-0) integrated intensity 
map of our observation field for M33. It reveals the 
clumpy nature of the molecular gas distribution. The 
size of these clumps ranges from < 100 pc (compara- 
ble to our beam size) to 500 pc, which corresponds to 
that of GMCs and/or GMAs. The northern arm that 
can be traced from (a, 5) = (li^33™40% -f30°46'0") to 
(ll^34'"30^ +30°49'0") is prominent. 

The surface mass density of molecular hydrogen was 
estimated by applying the CO-to-H2 conversion factor, 
Xco = 3 X 10^0 cm-2 (K kms^i)-! ijWilson fc Scovilld 
[1991 as follows: 

VMopc-V ■ V3 X 1020cm-2(Kkms-i)-V 

where i is the inclination angle of the galaxy [i = 51°). 

3.2. Star Formation Rate 

We calculated the SFR using the relation between the 
extinction-c orrected Ha line emi ssion and the SFR, as 
presented in lCalzetti et al.l (j2007l ): 

SFR(M0 yr-i) = 5.3 x lO'^^ [L{HaUs (2) 
-h(0.031±0.006)L(24^m)], 

in turn leading to 

SsFR(Moyr-ipc-2)=5.0 x 10"^ [I](Ha)obs (3) 
+ (0.031 ± 0.006)E(24 /im)] cos z. 



Breakdown of K-S law at GMC scales in M33 



3 



where the luminosities are in ergs ^ and L(24/im) is 
expressed as v L{v). The unit of I](Ha)obs and S(24 /im) 
is ergs~^ cm~^. 

We use the Ha image of M33 as given by 
iHoopes fc WalterbosI (12000) . It was obtained using the 
0.6-m Burrell- Schmidt telescope at the Kitt Peak Na- 
tional Observatory. The dimensions of the CCD are 
2048 X 2048 with a pixel size of 2". 03 and the total field 
of view is approximately 70' x 7 0'. The sensitivity is 
0.8 X 10-"ergs-icm-2arcsec-2 (jHoopes et al.ll200'l[ ). 

We retrie ved the 2 4- /im Multiband Imaging Photome- 
ter (MIPS: iRieke et al. 2004) datasets (AORs 3647744, 
3648256, 3648512, 3648768) of the basic calibrated data 
(BCD) created by the Spitzer Science Center (SSC) 
pipeline from the Spitzer Space Telescope (jWerner et al.l 
[2004) data archive. The mosaics were assembled by gath- 
ering all the BCDs except the first data frames of each 
observation because they have shorter integration times. 
The indiv idual calibrated frames were processed using 
MOPEX (|Makovoz fc Khanl[2005h for cosmic-ray rejec- 
tion and background matching was applied to overlap- 
ping fields of view. The final dimension of the each 
mosaic was 21' x 89'. After we removed the zodiacal 
light, and the noisy pixels near the edges by trimming, 
the four images were aligned with their coordinates and 
combined. The FWHM of point spread function is 5". 7. 

The Ha and 24-/im images were convolved into the 
same angular resolution as the C0( J= 1-0) data (19". 3), 
regridded to 7". 5 per pixel, and trimmed to match the 
CO ( J= 1-0) map. The rms noise is 5 x 10~^ erg s""'^ cm~^ 
in the Ha map, and 3 x 10^''' ergs^^ cm^^ in the 24- 
/xm map. The resultant SFR image is shown in Fig. 2 
with the overlaid CO ( J= 1-0) contours. The rms noises 
in the Ha and 24-/im maps result in a Esfr error and a 
sensitivity limit of 5 x 10~^^ Mq yr~^pc~^. 

Fig. 2 reveals a striking variety of star formation prop- 
erties. In regions where CO emission is detected, there is 
a large dispersion in the values of Ssfr of over four or- 
ders of magnitude. For example, the most massive GMC 
at (a, 5) = (l'^33'°9^6,-H30°49'7".3) has little star for- 
mation activity despite the large amount of molecular 
gas it contains. In the three major star-forming regions 
(NGC 604, NGC 595, and IC 133), there is a considerable 
difference in the amount of molecular gas associated with 
these star-forming regions. NGC 604 has its associated 
GMCs at the same position, NGC 595 has the ones offset 
from it, and moreover, none of the GMCs are associated 
with IC 133. 

4. DISCUSSION 

Fig. 3 presents the relationship between the SFR and 
molecular hydrogen surface densities for four angular res- 
olutions: 19".3 (-80pc), 60" (-240pc), 120" (-500pc), 
and 240" (~ 1 kpc). The panel of ^ 80 pc resolution is 
the K-S law plot with the highest spatial resolution for 
an extra-galaxy to date. Fig. 3 shows that apparent cor- 
relations exist between a-nd Ssfr at the ~ 1 kpc 
resolution, as already found in the disk-averaged data of 
M33 (Heyer et al. 2004). The best least-squares fit for 
the ~ 1 kpc resolution data with > 2(7 is 

log SsFR = (1.18 ± 0.11) log - (9.38 ± 0.05). (4) 

We have fitted the power law only to the data at the 1- 
kpc resolution, because the data at the other resolutions 



have a significant selection bias due to the omission of 
data points having < 2cr. The selection bias makes 
the correlation slopes appear steeper than real. 

Fig. 3 also shows that the correlation becomes looser 
with higher resolution, and it is hardly visible in the plot 
for ^ 80 pc resolution. Although most values lie 
within 10-40 Mq pc~^ where > 2(t, Ssfr values 
exhibit a wide range of over four orders of magnitude, 
from < 10~^° to ~ 10~^ M© yr"^ pc"^, indicating that 
the K-S law becomes invalid at this resolution. That is, 
the K-S law is valid only for averaging SFR and gas mass 
in large scales of several hundred parsecs. Our resolution, 
~80 pc is smaller than these scales. 

We have examined the effects of stochasticity on the 
estimatiion of SFR due to small sampling at high spatial 
resolutions. The initial mass function is not fully popu- 
lated when smaller regions that contain only a few stars 
in clusters are sampled. Thus, in regions with weaker 
extinction-corrected Ha emission, this effect may lead 
to significant scatter in the estimated Ssfr- We have 
estimated the significance of this effect on the 80-pc res- 
olution data. For instance, Ssfr = 10^^ Mq yr~^ pc~^ 
corresponds to the total flux of extinction-corrected Ha 
emission having i(Ha)corr — 6.6 x lO'^^ergs"^ within 
the 80-pc beam. Assuming a gas temperature of 10^ K 
and an electron density of 100 cm~^ in Hii regions, 
the ionizing photon flux is given by Q(H°) = 7.3 x 
10"L(Ha)corrS-i (IKennicuttlll988t IBrocklehurstl[T97l . 
Therefore, the corresponding total ionizing photon flux 
is Q(H°) = 4 9 X lO'^^s"^ According to Fig. 3 in 
iCervino et al.l (|2002f ). (5(H°) per unit mass varies with 
cluster age from 10**^ to 10"*^ s~^ ^o^- The ratio of the 
error in Q(H°) to the mean value, cr[(3(H*')]/Q(H°), is 
~ 10M"i/2 (the ratio can vary by a factor ^ 2 de- 
pending on the cluster age), where M is the cluster 
mass. Thus, the cluster mass required to produce the 
ionizing flux mentioned above is estimated to be 50 to 
5 X 10" Mq, which results in a[Q{Y{°)]/Q{Y{^) « 1 to 
0.04. In summary, in the 80-pc scale, the error due to 
the effect of stochasticity at Ssfr = 10~® Mq yr~^ pc~^ 
is K, 0.04-1 X 10~^ Mq yr~^ pc~^. The error decreases 
with larger Esfr and lower spatial resolutions. For ex- 
ample, the maximam cr[(5(H")]/Q(H°) in the 80-pc scale 
is 1, 0.4, and 0.1, and in the 240-pc scale it is 0.5, 0.1, 
and 0.05 for Ssfr = 10"^ 10-^ and IQ-^Mgyr-ipc-^, 
respectively. Although the error is signiflcant at smaller 
SsFR at the 80-pc resolution, it cannot explain the scat- 
ter of over four orders of magnitude in Ssfr- 

Now, we consider the possible causes for the breakdown 
of the K-S law at high spatial resolutions. In grand design 
spiral galaxies, we often see a systematic offset between 
the molecular gas arms and the star forming arms as a 
result of the density wave and the time delay between the 
accu mulation of gas an d its ionization by newborn stars 
(e.g.. lEgusa et al.ll2009[ ). If such an offset exists in M33, 
it may result in the breakdown of the K-S law when we 
observe the galaxy with a spatial resolution comparable 
to the offset. However, there is no such systematic offset 
between the CO and the star-forming arm in M33, as 
shown in Fig. 2. Therefore, the offset can be ruled out 
as a cause. 

In the Milky Way, young clusters recede from their 
parent GMCs with velocities of at least ~ 10 kms~^ 
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(|Leisawitz et al.|[l989[ ). As they evolve, the separation 
between clusters and their parent GMCs becomes signif- 
icant and is expected to be ^ 100 pc at a cluster age of 
^ 10 Myr. This drift scale is comparable to our resolu- 
tion, ^ 80 pc, and thus, the separation could well lead 
to the scatter in our plot. However, the existence of Hii 
regions without neighboring GMCs and GMCs without 
neighboring Hii regions in M33, such as the ones men- 
tioned in the previous section, cannot be explained by 
drift alone. 

The other cause is the difference in evolutionary stages. 
As we have shown in Fig. 2, there is a large diversity in 
the star formation activity among GMCs. Some GMCs 
in the Milky Way are at least p artially consumed an d 
dissociated by newborn clusters ([Leisawitz et al.|[l989D . 
Such GMCs are shown in Fig. 2 (i.e., the GMC at (a, S) 
= (l'^34'"17", 30'^52"00")). The diversity among GMCs 
arises from the evolutionary stages that each GMC goes 
through, as has been s hown in the mos t active star 
forming region NGC 604 ()Miura et al.|[2010l) . iChen et al.l 
([2010) studied the relation between Esfr and Ehi-i-H2 
in LMC for six GMCs in different stages of evolution, 
which are determined from the association of the under- 
lying stellar population. They also found a large scatter 
in SsFR, which exhibits a range of two orders of mag- 
nitude whereas the Eh i+h^ lies within a factor of two. 
iKawamura et al.l ()2009D classified molecular clouds in the 
LMC into three types according to the activities of mas- 
sive star formation: those with no signature of massive 
star formation, those with relatively small Hii region(s), 
and those with both Hii region(s) and young stellar clus- 
ter(s). They concluded that there is a sequential relation 
between these three types. That is, molecular clouds 
are first formed; then, stars are formed in the molecu- 
lar clouds; and finally, the newborn stars consume and 
dissociate the surrounding gas. This scenario is also ap- 
plicable to M33, as can be shown in NGC 604 by compar- 
ing it with warm/d ense gas traced by CO (/= 3-2) lines 
(jTosaki et al.ll2007D . 

Note that the extinction-corrected Ha emission is not 
really an appropriate tracer of SFR at resolutions as high 
as ~ 80 pc, because there is a spatial displacement be- 
tween ioni zing stars and the surroun ding Hii gas. Ac- 
cording to lRelano fc Kennicutti ()2009[ ). the displacement 
is ~ 10 pc to 50 pc in M33, and there is a radial gradi- 



ent of emission from the center to the boundary of Hii 
regions; in the order of ultraviolet. Ha, 24 /zm, and CO 
emission. Because of the radial gradient, the separations 
between ionizing stars and CO are larger than that be- 
tween Ha and CO. Thus, the "real" SFR would produce 
an even larger scatter in the K-S law plot for the ^80-pc 
resolution. The contribution of dust heating from older 
stars to the 24-/im emission can also be responsible for 
the scatter, although its effect is minor. 

Owing to our spatial resolution, which is comparable 
in scale to that of GMCs, we have clearly revealed the 
variety in the evolutionary stages among GMCs. The dif- 
ferent evolutionary stages of GMCs are the main cause 
of deviations from the K-S law. The higher the spatial 
resolution, the more individual star-forming regions and 
molecular clouds we can resolve, and the larger the de- 
viations seen in the K-S law plots. 

5. SUMMARY 

We conducted on-the-fly mapping of the northern part 
of M33 (30' X 20') in ^^CO ( J= 1-0) lines with the 45-m 
telescope at the Nobeyama Radio Observatory (NRO). 
We tested the relationship between the surface density 
of molecular gas mass and the SFR in an external galaxy 
with the highest resolution 80 pc) to date. We found 
that the star-forming activities in the GMC scale exhibit 
a wide range of over four orders of magnitude at positions 
where CO is significantly detected, such that Ssfr ^ 
10^^° to ~lO~^M0yr~^ pc~^, whereas their values 
are mostly in the range of 10-40 Mqpc"^. We also found 
that the K-S law of H2 gas, Esfr oc S^^, becomes invahd 
at GMC scales (-- 80 pc). The breakdown of the K-S 
law is attributed to the difference in the evolutionary 
stages of GMCs and the drift of newborn clusters from 
their parent GMCs that become apparent at this spatial 
resolution. 
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Right Ascension (J20CI0) 



Fig. 1. — (Left) BVHa image of M33 taken with SUBARU Suprime-Cam: courtesy of V. Vansevicius, S. Okamoto, and N. Arimoto. 
The rectangle represents our mapping area of CO(J=l— 0) (30' X 20'). (Right) Smoothed and clipped integrated intensity map of the 
CO (J= 1—0) lines. The contour levels are 1, 2, 4, 6, and 8 K kms~^, respectively. 
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Fig. 2. — CO integrated intensity map (contour) overlaid on the SFR map convolved into 19". 3 resolution. The contour levels are 1, 2, 
4, 8, and 16 K kms~^. The circle at bottom left represents the beam size. 



Breakdown of K-S law at GMC scales in M33 
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Fig. 3. — Star formation rate per unit area (Ssfr) versus surface density of H2 gas (Sh?) for four different resolutions: (a) 19". 3 (~80pc), 
(b) 60" (~240pc), (c) 120" (~500pc), and (d) 240" (~1 kpc). The broken red lines in each panel represent T.H2 = 2o- of the maps. The line 
in the lower-right panel represents the best least-squares fit to the ~ 1 kpc resolution data; log EsFR = (1.18 ±0.11) log Sjjj — (9.38 ± 0.05). 
The data points with < 2(7 are not used for the fitting. 



